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■<^j- , Abstract 

Uncertainties in predicting the nonlinear clustering of matter are among the most 
serious theoretical systematics facing the upcoming wide-field weak gravitational 
lensing surveys. We estimate the accuracy with which the matter power spectrum 
will need to be calibrated in order not to contribute appreciably to the error budget 

Q-f for future weak lensing surveys. We consider the random statistical errors and the 

systematic biases in P(k), as well as some estimates based on current N-body sim- 

£3 . ulations. While the power spectrum on relevant scales (0.1 < /c/ZiMpc -1 < 10) is 

currently calibrated with N-body simulations to about 5-10%, in the future it will 
have to be calibrated to about 1-2% accuracy, depending on the specifications of the 
survey. Encouragingly, we find that even the worst-case error that mimics the effect 
of cosmological parameters needs to be calibrated to no better than about 0.5-1%. 
These goals require a suite of high resolution N-body simulations on a relatively fine 
grid in cosmological parameter space, and should be achievable in the near future. 
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1 Introduction 



Weak gravitational lensing (WL) of distant galaxies is determined by the to- 
tal matter overdensity along the line of sight. In more technical terms, weak 
lensing measures the matter power spectrum weighted with a function that 
depends on the geometry of the lens-deflector-observer system. Therefore, in 
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order to predict the WL signal for a given set of cosmological parameters, one 
needs to be able to accurately calibrate the matter power spectrum P(k). 

The sensitivity of weak lensing surveys peaks at scales of about 10 arcminutes 
on the sky, corresponding to intervening structures of size ~ lh^ 1 Mpc at 
redshift ~ 0.5, which is approximately where the lensing efficiency of distant 
galaxies peaks (e.g. Jain & Seljak 1997). Therefore, the scales probed are 
nonlinear and need to be calibrated by N-body simulations. Unfortunately, 
one cannot afford to restrict the measurements to linear scales for which there 
are accurate analytical predictions (and which corresponds to scales > 1 deg. 
at z ~ 0.5), since the information from large scales only is poorly constraining 
due to cosmic variance (e.g. Huterer 2002). 

Therefore, accurate calibration of the fully nonlinear matter power spectrum 
P(k) is necessary to ensure that the cosmological parameter determination 
from weak lensing surveys is not compromised. Current accuracy in the matter 
power spectrum at mildly nonlinear scales is about 5-10% (White & Vale 
2004; Heitmann et al. 2004) which is sufficiently small for current surveys (for 
reviews of weak lensing, see e.g., Bartelmann & Schneider 2001, Van Waerbeke 
& Mellier 2003, Refregier 2003). However, the accuracy will need to be better 
for future wide-field surveys. The accuracy and agreement of various N-body 
simulations is rapidly improving and the currently popular fitting formulae for 
the power spectrum (e.g. Ma et al. 1999, Jenkins et al. 2001, Smith et al. 2003) 
will soon be insufficiently accurate and flexible and one will need to resort to a 
full suite of N-body simulations, interpolating on a grid of cosmological models 
(and thereafter, perhaps, building a new set of more flexible and accurate 
fitting formulae). 

In this work we consider the requirements of future deep, wide-field WL sur- 
veys such as Dark Energy Survey 1 , PanSTARRS 2 , Supernova/Acceleration 
Probe (SNAP 3 ; Aldering et al. 2004) and Large Synoptic Survey Telescope 
(LSST 4 ). For the fiducial model we assume sky coverage of 1000 sq. deg. with 
100 galaxies per arcmin 2 . We vary these numbers later to determine the ro- 
bustness of our results to survey specifications. We assume a flat universe with 
matter energy density relative to critical Qm = 0.3, dark energy equation of 
state w — — 1, and power spectrum normalization cr 8 = 0.9. We use the spec- 
tral index and physical matter and baryon energy densities with mean values 
n = 1.0, VL M h 2 = 0.140 and Vt B h 2 = 0.023 respectively. The parameters whose 
future accuracy can most benefit from future WL surveys are Qm, w and as, 
and we pay particular attention to these three. We do not add any priors 
whatsoever, since we are considering the efficacy of WL surveys in their own 

1 http://cosmology.astro.uiuc.edu/DES 

2 http://pan-starrs.ifa.hawaii.edu 

3 http:/ / snap.lbl.gov 

4 http://www.lsst.org 
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right; moreover, we have checked that reasonable priors (for example, those 
based on expected accuracies from the Planck experiment) have minimal effect 
on the three aforementioned parameters and all of our results. Note that we 
do not consider the time variation in w since the best-measured mode of an 
arbitrary w(z) is about as well measured asw = constant (Huterer & Stark- 
man 2003), and it is this particular mode, being the most sensitive to theory 
systematics, that will drive the accuracy requirements. Hence, it is sufficient 
to consider the constant w case. Throughout we consider lensing tomography 
with 10 redshift bins equally spaced between z = and z = 3 and use the 
lensing power spectra on scales 50 < £ < 3000. Over this range of multipoles, 
statistical properties of the lensing fields are nearly Gaussian and, furthermore, 
complex baryonic effects that are expected to be significant on smaller scales 
are strongly suppressed (e.g. White 2004, Zhan & Knox 2004). [Cosmological 
parameter extraction from WL tomography has been extensively discussed in 
the past, e.g. Hu 1999, Huterer 2002, Hu 2003, Refregier et al. 2003, Takada 
& Jain 2004, Takada & White 2004, Song & Knox 2004, Ishak et al. 2004.] 

The convergence power spectrum at a fixed multipole £ and for the ith and 
jth tomographic bin is given by 



™ iz WM P ( ' A, (!) 
r(z) z H(z) \r{z) J 



where r(z) is the comoving angular diameter distance and H(z) is the Hubble 
parameter. The weights Wi are given by Wi(x) — § Hq gi(x) (1 + z) where 
9i(x) = r (x) J™ dx s ni(Xs)r(Xs ~ x)/r(Xs), X is the comoving radial distance 
and rii is the comoving density of galaxies if Xs falls in the distance range 
bounded by the ith redshift bin and zero otherwise. We employ the redshift 
distribution of galaxies of the form n(z) oc z 2 exp(— z/zq) that peaks at 2zo = 
1.0 (for LSST, we use 2z = 0.7). The observed convergence power spectrum 
is 

cm= p m+^ {J ^ L ^ ( 2 ) 

where {'jf nt ) 1 ^ 2 is the rms intrinsic shear in each component which we assume 
to be equal to 0.22, and fii is the average number of galaxies in the ith redshift 
bin per steradian. The cosmological constraints can then be computed from 
the Fisher matrix 

F^E^g, (3, 

where Cov 1 is the inverse of the covariance matrix between the observed 
power spectra whose elements are given by 



c*w, cm] = {2e+1 " f M [cmcm + cmcm] ■ w 



Cov 

assuming that the convergence field is Gaussian, which is a good approxima- 
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Parameter 


SNAP Error 


LSST Error 


^DE 


0.008 


0.003 


w 


0.052 


0.025 


0"8 


0.006 


0.003 



Table 1 

Fiducial 1-a errors in Q.m, w and ag for the SNAP-type and LSST- type surveys 
assuming no systematics. These errors were obtained by marginalizing over all (five) 
other cosmological parameters. 

tion at £ < 3000. We alternatively consider a SNAP-type survey (1000 sq. deg., 
100 gal/arcmin 2 ) and an LSST-type survey (15000 sq. deg., 30 gal/arcmin 2 ). 
The fiducial marginalized accuracies in Qm, w and erg for the two surveys, 
without any theoretical systematics, are given in Table 1. 

In Section 2 we assume random errors in the matter power spectrum P(k) and 
determine their maximal allowed size so that they do not appreciably add to 
the overall error budget in future surveys. In Section 3 we explore biases due 
to several specific sources of error, and also consider the worst-case scenario 
when errors in P{k) conspire to maximally bias the cosmological parameters. 
We conclude in Section 4. 



2 Statistical errors in P(k) 

First, let us parameterize the deviations in the matter power spectrum around 
its fiducial value by 30 parameters 5 In Pi in bins spaced equally in log k from 
k = 0.1/iMpc -1 to k — lO/iMpc -1 . More precisely, if k is in the ith bin, we 
allow the variations 

P(k,z) P(k,z){l + 5\nP i ) (5) 

where 5 In Pj are assumed to be independent gaussian random variables, not 
determined by the cosmological parameters, with mean zero. Note too that 
we assume that relative variations in P(k,z), parametrized by £lnPj, are 
independent of z. This simplifying assumption is motivated by the fact that the 
source of error in In P(k, z) is imperfect accuracy of N-body simulations, where 
biases in In P(k, z) are likely to be weakly dependent on z in the range relevant 
for WL. We compute the linear power spectrum using the fitting formulae 
of Eisenstein & Hu (1999) which were fit to the numerical data produced by 
CMBFAST (Seljak & Zaldarriaga 1996). We generalize the formulae tow ^ -1 
by appropriately modifying the growth function of density perturbations. To 
complete the calculation of the full nonlinear power spectrum we use the fitting 
formulae of Smith et al. (2003). 
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Fig. 1. Effects of random statistical error on the power spectrum P(k). In the top 
panel we represent P(k) in 30 bins distributed uniformly in log A; from k = 0.1 
to 10/iMpc -1 and show the unmarginalized accuracy in each band power, 5lnP{, 
assuming that the cosmological parameters are perfectly known. It is clear that the 
scale that is most sensitive to weak lensing measurements is k ~ 1/iMpc -1 , and 
therefore P(k) around this scale needs to be calibrated most accurately. We also 
find that about 10 band powers of P(k) in the plotted range are independent. The 
bottom panel shows the increase in the error in cosmological parameters TO , w and 
eg as a function of prior imposed on each of the 10 nuisance parameters 5 In Pi . 
The solid and dashed curves are the results for SNAP- and LSST-type surveys, 
respectively. The required accuracy is about 2% from erg and 5-10% from &m and 
w for SNAP, and about a factor of two more stringent for LSST. 
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We assume for a moment that all cosmological parameters are precisely known 
and fixed. The top panel of Figure 1 shows that the unmarginalized accura- 
cies in 5 In P; L range from a few hundredths to a few tens, except on very small 
scales. Not surprisingly, the best accuracy is achieved on scales of ~ l/iMpc -1 
to which the lensing power spectra of our interest are most sensitive. From 
the covariances between different 5 In Pi (i.e. the off-diagonal components of 
the Fisher information matrix), we find that only about 10 bins in P(k) are 
independent - in other words, the weak lensing kernel is broad enough to 
differentiate only around 10 parameters in P(k) in this range of k. We then 
investigate how accurately these 10 numbers need to be known for the cos- 
mological parameter accuracies not to be affected appreciably. We add equal 
priors to 10 values <51nPj (i = 1 . . . 10) in bins spaced equally in log A; from 
k = 0.1/iMpc" 1 to k — 10/iMpc" 1 and marginalize over them, then com- 
pute the degradation in the marginalized accuracies of Qm, °"s an d w. The 
bottom panel of Figure 1 shows that a prior on the P(k) values of a few per- 
cent (2% for as, 10% for Qm, 5% for w) would render the power spectrum 
essentially "precisely known" for a SNAP-type experiment (1000 sq. deg., 100 
gal/arcmin 2 ), and enable determination of the three cosmological parameters 
without degradations greater than about 10-20%. The required accuracies are 
about two times more stringent for an LSST-type experiment (15000 sq. deg., 
30 gal/arcmin 2 ). Also interesting to note is that as is more degraded than Q m 
and w. This can be easily understood since, while a$ depends on the lensing 
power spectra solely through its dependence on P(k), the parameters Q m and 
w also enter the lensing geometrical factor that only depends on the cosmic 
expansion history. In other words, most information on dark energy param- 
eters comes from the cosmic geometrical factor - provided that the redshift 
distribution of galaxies is precisely known (see Ma, Hu & Huterer 2004 and 
Huterer et al. 2004 for a discussion of the required accuracy of photometric 
redshifts for weak lensing tomography). 



3 Systematic biases in P(k) 



We now estimate the bias in cosmological parameters, Sp i} due to the system- 
atic bias in P(k). We can easily estimate Spi using the Fisher matrix formalism. 
Let Pi be the true values of cosmological parameters, and C£ and be the 
true and perturbed values of the convergence power spectra that result from 
the bias in P(k). For clarity we have denoted a pair of redshift bins by a sin- 
gle subscript a. Then, as long as the perturbations in the convergence power 
spectra are small compared to their measurement accuracy, one can derive the 
expression to estimate the bias in the cosmological parameters: 



c K a ( 



Cov" 



c K a ( 



(6) 
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where F^ 1 is the inverse of the Fisher matrix for the parameters pi, and the 
summation over j, a and (5 is implied. 

3.1 Some examples of realistic biases in P(k) 

We first estimate the bias in cosmological parameters due to the systematic 
error in P{k) estimated in N-body simulations. As mentioned earlier, the only 
reliable way to make accurate theoretical predictions for WL is to rely on nu- 
merical simulations and perform ray-tracing photons through simulated large- 
scale structures (e.g., Jain, Seljak & White 2000, White & Hu 2000, Heitmann 
et al. 2004). The accurate model predictions for a given cosmological model 
then have to be constructed from a sufficient number of the simulated lens- 
ing map realizations to reduce the statistical variance. However, clearly the 
simulated lensing map will contain uncertainties due to various numerical 
limitations. Vale & White (2003) derived a useful, approximate expression for 
how the finite resolution of the simulation and shot noise due to finite particle 
number affect the 3D matter power spectrum: 

A»(Jb) = (a^K^ 2 + e-^ 2 , (7) 

where Ag(fc) = k 3 P e (k) / (2n 2 ) is an effective 3D matter power spectrum, 
A-m(k) is the underlying true power spectrum, n is the mean particle num- 
ber density and a n and a g are characteristic resolution limits of the N-body 
simulation and the Fourier grid of the ray-tracing simulation, respectively. Fol- 
lowing Vale & White (2003) we assume that a g and o~ n can be approximated 
as a g = 0.54Lbox/-/Vgrid and a n = 0.05n -1 / 3 , where iVg rid is the ray-tracing grid 
number and L^, ox (Mpc) is the box size of N-body simulation used. Note that 
n is expressed in terms of L box and N-body particle number N as n = N/L 3 ox . 

The top panel of Figure 2 shows the fractional differences between the conver- 
gence power spectra with and without the model numerical error in P(k). 
We assumed simulation parameters of Lbox = 200Mpc//i, N = 512 3 and 
-Vgrid = 4096, and in this case the term including a n in Eq. (7) which models 
the finite force resolution in the N-body simulation provides dominant con- 
tribution to the error in P(k) (also see Figure 3 in White & Vale 2004). We 
find that the model error in P(k) causes the suppression of the lensing power 
spectrum amplitudes by £ 10% for £ < 3000. Since a given angular mode £ 
corresponds to smaller k as redshift increases, the lensing spectra of higher 
redshift bins arise from density perturbations of smaller k and are therefore 
less affected by the numerical effects. The shaded region shows the la statis- 
tical error in the lensing spectrum measurement for the fiducial survey, where 
the error arises mainly from the cosmic sample variance at t £ 10 3 , while the 
shot noise due to the intrinsic ellipticities is dominant at larger £. The depar- 
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Fig. 2. The top panel shows fractional differences between the lensing power spectra 
with and without the model biases in P(k) due to finite size and resolution of the 
simulation, a la Eq. (7). The solid curves show the auto spectra of 10 redshift 
bin tomography, where the redshift bins are equally spaced between z = and 
z = 3 and we assumed the simulation parameters of L = 200 Mpc//i, n = 512 3 , and 
-^grid = 4096. The shaded region shows the la statistical error for measuring the 5th 
redshift bin spectrum, C55, for the SNAP-type survey (1000 deg 2 , 100 gal/arcmin 2 
and (7int) = 0.22). Bandwidth of Al/l ~ 0.11 is assumed. The bottom panel shows 
how the bias in P(k) propagates into the bias in cosmological parameters. We plot 
the bias in the cosmological parameters relative to their marginalized error vs. the 
ray-tracing simulation grid number N^^. The solid and dashed curves show the 
results for the two cases of N-body simulation particle number, ./V = 256 3 and 
iV = 512 3 , while the N-body simulation box size is fixed to 200 Mpc/h. 
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tures in the lensing spectra of low redshift bins are larger than the statistical 
errors even though we only consider the measurements at t < 3000. 

Thus, it is important to investigate how the numerical inaccuracy in P(k) 
propagates into the biases in cosmological parameters, and we address this 
in the bottom panel of Figure 2. As with the statistical error in P(k) (see 
Figure 1), the biases in P(k) lead to a larger bias in og than those in Q m 
and w. It is reassuring to see that the simulations with iVg ri d > 4096 and 
N = 512 3 , which are feasible with current numerical resources, lead to the 
bias in w and Q m of less than about 20% relative to the marginalized error. 
Moreover, in the future one can expect to be able to calibrate, and partly 
correct for, the remaining bias in simulations by using a small number of even 
higher resolution simulations as a reference. 

The effect of baryons on the non-linear mass clustering is another potentially 
harmful systematic effect on P(k) (White 2004, Zhan & Knox 2004) and its 
calibration requires hydrodynamical N-body simulations that include baryon 
dissipation. Baryons, which constitute more than 10% of the total mass, cool 
and condense within a halo, and their effects are most pronounced for halos 
of mass £ 1O 13 M (White 2004). The X-ray observations indicate that the 
hot baryons are likely to have a smoother profile than dark matter. Zhan & 
Knox (2004) estimated the effect of hot baryons on the lensing power spec- 
trum assuming an isothermal /3-model to describe the density profile of the 
baryon within the host halo. Their fiducial model leads to modifications in the 
lensing power spectrum of a few percent around £ ^ 10 3 for source redshift 
z = 1 (see Figure 2 in their paper). This effect is of similar amplitude as that 
shown in Figure 2, suggesting that the presence of hot baryons does not neces- 
sarily lead to drastic biases in the cosmological parameters. However, a much 
more extensive investigation that uses a suite of hydrodynamical simulations 
is needed to properly assess the effect of baryons on cosmological parameter 
estimation. 



3.2 Worst-case Systematic Bias in P(k) 

The most damaging systematic error in calibrating P(k) will generally be error 
that mimics the behavior of cosmological parameters. Somewhat confusingly, 
however, the worst-case biased P(k) (call it Q(k)) is typically not the matter 
power spectrum for some neighboring values of cosmological parameters. For, 
Q(k) is weighted with the geometric factor to produce the observable conver- 
gence power spectrum, and the geometric factor is known exactly (provided 
that the redshift distribution of galaxies is known). Instead, the worst-case 
Q(k) will be such that, when weighted by the geometric factor of the true cos- 
mological model, it produces the convergence power spectrum which precisely 
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matches some neighboring cosmological model. 



While finding the actual worst-case Q{k) is a thorny problem, we can compute 
the desired results in the approximation that the perturbation in P(k) is small. 
The assumed deviation in the convergence power spectra is solely due to the 
bias in the matter power spectrum 

N 

C K a (l)-C K a (£) = Y,W(k a )5P(k a ), (8) 

a=l 

where we have discretized the integral as the sum over N lens planes, or 
equivalently over N values of the wavenumber k = £/r(z), and defined W(k) 
appropriately; c.f. Eq. (1). Then Eq. (6) can be rewritten as 

N 

6p i = 'E(6]nP(k a ))c®, (9) 

o=l 

where we defined 

dC K (£) 

ch^E^y^WCov" 1 [CZ(£),C$(£)] (10) 

where the summation over j, a and /3 is implied. Let us assume for a moment 
that the bias in lnP(fc) is known to be smaller, by absolute value, than some 
fixed number Sin P. Then Eq. (9) shows that the worst-case bias in P(k) - 
one that maximizes Spi - takes the value ±(51nP) on each lens plane, with 
the sign equal to the sign of on that plane. The resulting worst-case bias 
in the cosmological parameter is Spi = (<51nP) Yla=i \ c $\- 

Turning the argument around, if the bias in a given cosmological parameter 
Pi is to be smaller than Sp^ it is sufficient to trust the matter power spectrum 
to an accuracy better than 

S\ n P = (11) 

and this is true irrespective of the conspiring shape of 5\nP(k). Figure 3 
shows the required accuracy in P(k), 5\nP, assuming the worst-case scenario 
where the bias equals the la accuracy in pj. The required accuracy is shown as 
a function of area of the survey for two values of the density of galaxies (100 
and 30 galaxies/arcmin 2 ) and for three cosmological parameters that are of 
greatest interest to us (£Im, °8 an d w). As before, we find that more ambitious 
surveys have more stringent requirements on the systematic control of P{k). 
We find that, in order for the worst-case systematics to be subdominant to 
the la error bars in the three cosmological parameters, a survey like SNAP 
will need a 1% calibration in P(k) over the relevant scales, while LSST needs 
a better than 0.5% calibration. We emphasize that these are the required 
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Fig. 3. The necessary accuracy in P(k) for the worst-case scenario of bias in P(k), 
defined so that it produces the maximal bias in the three cosmological parameters 
shown. The required accuracy, plotted as a function of the sky coverage of the 
survey, is defined so that it limits the biases in the cosmological parameters to be 
no larger than their la fiducial accuracies. We show the cases of n g = 100 and 30 
gal / arcmin 2 

sensitivities for the worst-case scenario of bias in P(k), where the bias in 
P(k) conspires to introduce the maximal bias in the cosmological parameters. 
Section 2 shows that a more common error in P(k) that is uncorrelated with 
cosmological parameters produces biases that are a factor of a few smaller. 



4 Conclusions 



Accurate knowledge of the full nonlinear matter power spectrum is crucial for 
the future weak gravitational lensing surveys to achieve their full potential. In 
this paper we have explored the required accuracy in P(k, z) so that it does 
not bias the cosmological parameter extraction from future wide-field surveys. 
This analysis complements other analyses of theoretical and observational sys- 
tematics in weak lensing measurements (e.g. Vale et al. 2004, Huterer et al. 
2004). 

We find that the power spectrum sensitivity is greatest at two decades in 



11 



wavenumber centered around k = 1 Mpc/h comoving. Calibration of P(k) 
on those scales needs to be better than a few percent (« 2% for cr 8 and 5- 
10% for fl M and w) for the SNAP-type experiment, and about a factor of two 
more stringent for the LSST-type experiment. We also consider the systematic 
biases in the power spectrum in N-body simulations. While the simulations' 
accuracy is already nearly sufficient, the cooling and collapse of baryons, and 
even the clustering of massive neutrinos (Abazajian et al. 2004) may lead to 
large errors if they are not modeled accurately. There exist a number of options 
to protect against such effects (Huterer & White 2004), but exploring those is 
beyond the scope of this work. The most stringent accuracy requirement on 
P{k) is obtained for the worst-case scenario where the bias in P(k) conspires 
to produce the convergence power spectrum that is precisely reproduced by 
cosmological parameters with biased values. Fortunately, we find that even 
this worst-case scenario leads to tolerable biases in the cosmological parame- 
ters provided P(k) is calibrated to 0.5-1% accuracy over the relevant scales. 
Note too that, in order to achieve ~ 1% accuracy in P(k), current analytic 
fits for the linear power spectrum are not accurate enough either, and one 
instead needs to numerically solve the coupled Einstein, fluid and Boltzmann 
equations as done e.g. in the CMBFAST code (Seljak & Zaldarriaga 1996). 

As we are entering the era of powerful wide-field WL surveys, currently popular 
fitting formulae for the power spectrum will soon become inadequate and one 
will need to resort to a full suite of N-body simulations, interpolating on a 
grid of cosmological models. First steps in such an approach were recently 
undertaken (White & Vale 2004) although many details remain to be worked 
out. We have shown here that the required accuracy is nearly within reach 
with current simulations, although the number of models that are necessary 
to be run might be very large. With the effort of the weak lensing community 
and ever-more powerful computers, the necessary simulation results should 
be in hand by the end of the decade, enabling weak lensing to realize its full 
power as a probe of dark energy and matter distribution in the universe. 

DH is supported by the NSF Astronomy and Astrophysics Postdoctoral Fel- 
lowship. We thank Kev Abazajian, Gary Bernstein, Salman Habib, Bhuvnesh 
Jain, Lloyd Knox, Eric Linder, Chris Vale and Martin White for useful con- 
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